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Abstract 

We have studied the magnetic properties of free cobalt clusters in two semi-independent Stern- 
Gerlach deflection experiments at temperatures between 60 and 307 K. We find that clusters 
consisting of 13 to 200 cobalt atoms exhibit behavior that is entirely consistent with superparam- 
agnetism, though complicated by finite-system fluctuations in cluster temperature. By fitting the 
data to the Langevin function, we report magnetic moments per atom for each cobalt cluster size 
and compare the results of our two measurements and all those performed previously. In addition 
to a gradual decrease in moment per atom with increasing size, there are oscillations that appear 
to be caused by geometrical shell structure. We discuss our observations in light of the two com- 
peting models for Langevin- like magnetization behavior in free clusters, superparamagnetism and 
adiabatic magnetization, and conclude that the evidence strongly supports the superparamagnetic 
model. 

PACS numbers: 36.40. Cg, 75.50.Cc, 75.50.Tt 



I. INTRODUCTION 



Interest in small magnetic systems predates the emerging field of nanoscience by several 
decades.-i^ Magnetism in small systems is effectively the frontier between atomic magnetism, 
surface magnetism, and bulk magnetism. Because of its strategic position, small system 
magnetism has both theoretical importance in elucidating how structural and electronic 
order evolve from one dimensional extreme to the other and practical importance in such 
technologies as magnetic data storage and high-performance magnetic materials. 

Magnetism in small and low-dimensional systems has been studied since the 1950s using 
powders, granular metals, bulk surfaces, and supported films and particles.-i^ Phenomena 
observed in these contexts included superparamagnetism and surface-enhanced magnetism. 
Superparamagnetism occurs when the overall moment of a small magnetically ordered par- 
ticle becomes thermally decoupled from its lattice so that it can respond paramagnetically 
to an applied magnetic field.-'^'^i^ Surface-enhanced magnetism results when the decreased 
coordination number and increased density of states at a surface reduces the quenching of 
spin and orbital magnetism relative to what occurs in the bulk.-i^ 

More recently, however, the magnetic properties of atomic clusters have been studied in 
the isolation of vacuum using molecular beam techniques.-'2iii2i>iiii^'i^i^'i^ii^'i^>i^ii2i>22^ That 
isolation has at least three interesting consequences: it frees the clusters from all chemical 
interactions at their surfaces, it decouples them from external heat baths, and it emphasizes 
their conserved quantities. As we will discuss, those isolation effects both simplify and 
complicate the observed behaviors of magnetic clusters in a beam. 

One might expect that a beam of magnetic molecules of spin S would, upon passing 
though a gradient magnetic field, split into 25+1 beamlets. This behavior is indeed the 
high-field case for atoms and small molecules.— One might also argue that the atomic 
moments /i in a ferromagnetic particle would couple together and orient themselves along the 
applied magnetic field, and that the particle as a whole would then accelerate and deflect 
toward high field in response to the force V(A^/^-B). That behavior is the high- field, 
low-temperature limit for single-domain condensed matter particles. The clusters we study 
are intermediate in size and we study them at moderate temperatures and magnetic fields. 
Anticipating their behaviors is therefore not so simple. 

The first magnetic deflection experiments on free cobalt clusters observed that these clus- 
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ters always deflect toward strong field by an amount that increases as the applied magnetic 
field increases, as the relevant temperature decreases, and as the number of atoms in the 
cluster increases.— In the present measurements, we again observe deflections that increase 
with increasing magnetic field and decreasing temperature, and that are entirely consistent 
with superparamagnetic behavior in which each cluster acts as its own heat bath at its vi- 
brational temperature.— The influence of cluster size, however, is more complicated than we 
observed originally. We now find that the each cobalt cluster's magnetic moment depends 
subtly on the number and arrangement of its atoms. 

In superparamagnetism, the effective magnetic moment per atom n^s is reduced from the 
true moment per atom n by the Langevin function C 



where N is the number of atoms in the particle, B is the applied magnetic field, ks is 
Boltzmann's constant, and T^ib is the vibrational temperature. This reduction occurs when 
thermal agitation decouples the particle's magnetic moment from its lattice and causes 
that moment to explore all possible orientations in rapid succession. Since the orientation 
fluctuations occur on a subnanosecond timescale,-'^'^ a millisecond-timescale measurement 
of the particle's magnetic momentum per atom will yield the time-averaged projection of 
l-i onto the applied magnetic field axis B — that is, the experiment will measure /iefr- This 
experimental moment per atom increases linearly with applied magnetic field for weak fields, 
but eventually saturates at /i. 

Equation ([T]) assumes that the total moment N/i is large enough to be treated classically, 
otherwise the Brillion function must replace the Langevin function. It also assumes that 
all orientations of the moment relative to the particle lattice are energetically equivalent, 
that the measurement time is long compared to the thermal fluctuation time so that the 
particle thoroughly explores the Boltzmann distribution during the measurement, and that 
the particle is coupled to an infinite heat bath at temperature Tvib.— 

While free clusters cannot meet these assumptions perfectly, there is considerable evi- 
dence that Eq. ([T]) is valid for cobalt and most other magnetic clusters over a considerable 
range of experimental conditions. We note here that an adiabatic mechanism of magnetiza- 
tion has recently been proposed.-^ That alternative theory arose as an explanation for the 
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broadening of the cluster beam as it deflects. Interestingly, the adiabatic mechanism report- 
edly gives the same Langevin-like reduction in measured magnetic moment that is predicted 
by superparamagnetic theory. We will discuss these two competing explanations later in 
this paper. For the present, we will simply interpret our experimental data using Eq. ([T]) in 
order to obtained values for fi, the magnetic moments per atom of the cobalt clusters, and 
make occasional reference to superparamagnetism, the model that we ultimately conclude 
is responsible for our experimental observations. 

What distinguishes the present study from the one performed more than a decade ago in 
our laboratory^^ is that we are now able to report /i for each cluster size individually. In the 
earlier study, our limited sample required us to average across cluster sizes. We therefore 
had to assume that /i was independent of cluster size, at least over the size range considered 
(40-200 atoms). In the present work, however, we have studied each cluster size individually, 
so that magnetic effects due to cluster size and structure can be recognized. 

The present study consists of two semi-independent experiments, which we will refer 
to as El and E2. These experiments were conducted before (El) and after (E2) a major 
renovation of the experimental apparatus. Most signiflcantly, they were performed with 
radically different cluster sources and data acquisition protocols. Although El and E2 share 
the same gradient-fleld magnet and mass-spectrometer, their semi-independence makes them 
a useful test of our ability to control cluster temperature and therefore make meaningful 
measurements. That the results of these two studies are so similar provides considerable 
support for our assertion that we understand the temperatures of clusters in our beams. 

II. EXPERIMENTAL PROCEDURE 

In concept, the experiment is essentially the Stern-Gerlach deflection technique applied 
to cobalt clusters. We produce a narrowly coUimated beam of temperature-controlled cobalt 
clusters, pass that beam through a gradient magnetic fleld, and measure the deflections and 
masses of the clusters at a distance downstream from the magnet. While we have described 
this technique previously^-, we discuss it again here primarily to point out the differences 
between the two experiments. El and E2. 

In both experiments, cobalt clusters were grown from atoms vaporized off a cobalt sam- 
ple by the focused second-harmonic light (532 nm) of a pulsed Nd:YAG laser. The laser- 
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produced plume of cobalt vapor was captured in a pulse of dense helium gas, where cluster 
growth and thermal equilibration occurred. This mixture of gas and clusters then became 
a cluster beam through a seeded supersonic expansion into vacuum. 

In experiment El, cobalt vapor and helium mixed in a cylindrical chamber approximately 
8 mm in diameter and 17.5 mm long. Clusters grew in this "waiting room" and came 
gradually into thermal equilibrium with its walls before expanding into vacuum through 
the chamber's conical exit nozzle. In practice, however, clusters sprayed out of the nozzle 
continuously once cobalt atoms were injected into the helium pulse. We found that only 
those clusters that remained in the source as long as possible and didn't emerge from the 
chamber until it was almost empty had time to reach thermal equilibrium with the chamber 
walls. The measurements we report for El were all conducted with clusters that had resided 
in the waiting room long enough (>1.5 ms) to reach equilibrium. The range of cluster 
temperatures used in El was 63 K to 307 K, with 63 K being approximately the coldest 
source temperature that the closed-cycle helium refrigerator could maintain. 

In experiment E2, cobalt vapor and helium mixed and cooled in a cylindrical channel 
2.5 mm in diameter and 150 mm long. This channel had approximately the same volume 
as the "waiting room" used in El, but a much larger surface area. Cluster growth occurred 
near the start of this channel and the resulting clusters were thermally equilibrated in the 
channel's final 100 mm, which is temperature-controlled. The channel ends with a 1 mm 
cylindrical nozzle, through which the mixture of vapors underwent a seeded supersonic 
expansion to form a cluster beam. The increased surface to volume ratio of this source, 
coupled with the fact that even the fastest clusters take more time to get from the point 
of vaporization to the nozzle, means that a greater fraction of the clusters arc in thermal 
equilibrium with the source upon exiting. All of the measurements we report for E2 were 
made using equilibrated cobalt clusters. The range of cluster temperatures used in E2 was 
60 K to 100 K, again hmited at the low end by the closed-cycle helium refrigerator. 

Both experiments shared the same beam shaping apparatus. The cluster beam travels 
through a skimmer and two narrow slits (0.4 mm wide by 2.5 mm high) that coUimate 
the beam, and through a chopper wheel that slices it longitudinally. The chopper wheel 
performs two important functions: it allows us to select only those clusters that have come 
into thermal equilibrium with the source and it permits us to measure the velocity of the 
clusters. The wheel rotates at 180 Hz and is open to the cluster beam for ~20 /is. After 
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passing through these colhrnating shts and chopper, the cluster beam is a tightly grouped 
packet, narrow in all three spatial dimensions. 

Clusters in this packet then pass through the gradient field magnet, where they accelerate 
in response to the force 'V{N^ • B). Since the magnet's field gradient is perpendicular to 
the cluster beam axis, magnetic clusters are defiected, and it is that defiection that we 
subsequently measure. 

The gradient field magnet is 250 mm long and its gradient is nearly constant in the region 
traversed by the cluster beam.— We can vary the field gradient experienced by the clusters 
from T/m to 360 T/m, while simultaneously varying the field at the position of the cluster 
packet from T to 0.951 T. 

After leaving the magnet, clusters pass through a ~1 m field- free drift region and are then 
ionized by a spatially filtered pulse of laser light (193 nm, ArF Excimer) in a time-of-fiight 
mass spectrometer. This narrow light beam propagates antiparallel to the cluster beam and 
ionizes only those clusters that have defiected a specific distance from the undefiected beam 
center. By recording the resulting mass spectra at a range of positions relative to the beam 
center, we obtain a complete defiection profile for each of the cluster sizes present in the 
beam. 

In experiment El, these defiection profiles were obtained one magnetic field at a time. 
With the source producing pulses of clusters at 20 Hz and the magnetic field and field 
gradient set to one value, we scanned the laser beam back and forth across the cluster beam 
and accumulated cluster deflection profiles for several hours. We then chose another field 
and field gradient and repeated this exercise. We periodically turned off the magnet to 
obtain a zero-field profile to establish the undefiected beam center. 

Because the experimental apparatus is ~3 m long and the temperature in the laboratory 
varies over the course of a day, there are small but unavoidable movements in the equipment. 
Furthermore, the intensity and direction of the cluster beam itself drifts with time as the 
sample ages and soot accumulates on the source's internal surfaces. These changes produce 
gradual shifts in the measured defiection profiles, and the one-field-at-a-time measurement 
approach used in experiment El makes it difficult to remove these shifts from the data. 

Conversely, in experiment E2, through computer automation, we collected data randomly 
in both position and magnetic field, with the magnetic field switching between various values 
(including zero) and the ionizing laser shifting to a new randomly selected position a few 
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times a minute. Although the source used in E2 produced pulses of clusters at 30 Hz, a 
single experimental run still spanned an entire day. Although the long duration of each 
experimental run in E2 may have increased the statistical uncertainty, experimental drift no 
longer contributed significantly to the systematic uncertainty. 

III. RESULTS 

A. Residence Time Studies 

Because clusters spend >100 /xs passing through the gradient field magnet, our experi- 
mental measurement timescale is long compared to the subnanosecond timescale of super- 
paramagnetic fluctuations. Assuming that the cobalt clusters behave superparamagnetically 
in these experiments, the magnetic moment responsible for a cluster's deflection by the gra- 
dient field magnet is the time-averaged projection of that cluster's true moment onto the 
magnetic field. The experimentally measured moment per atom /iexpt is thus the effective 
moment per atom /ieg given by Eq. ([1]). 

One of the best indications that cobalt clusters obey Eq. ([1]) (and that they behave su- 
perparamagnetically in our experiments) is that /^expt has the correct functional dependence 
on applied magnetic field B and cluster vibrational temperature T^ih- In order to be sure 
that we know T^ih, however, we must prove that the clusters have resided in the source 
long enough to reach thermal equilibrium with its walls. We can then use the source's wall 
temperature as T^ib in all the subsequent analysis. 

Because knowing Tvib is so important, we regularly perform studies in which we measure 
/iexpt as a function of residence time in the source. The residence time for a particular group 
of clusters is defined as the interval between the Nd:YAG laser pulse and the moment those 
clusters emerged from the nozzle into vacuum. We can determine the moment at which a 
particular group of clusters left the nozzle from the timings of the chopper wheel and the 
ionizing laser pulse and from the distances separating these components from the nozzle. 
Knowing when the clusters left the nozzle allows us to calculate their residence time to 
±10 /is. 

As shown in Fig. [1], we always observe an initial rise in /iexpt with increasing residence 
time, followed by a long period of saturation in which additional residence time has no effect 
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FIG. 1: The experimental magnetic moment per atom /iexpt initially increases with residence 
time in the source but eventually saturates when the cluster vibrational temperature r^ib reaches 
equilibrium with the source temperature, (a) is a saturation curve from experiment El, taken 
for Coioo at a source temperature of 123 K and a magnetic field of 0.269 T, and (b) is one from 
experiment E2, taken for Cogs at a source temperature of 100 K and a magnetic field of 0.204 T. 

on the measured moment per atom. The timescales for the waiting room source used in 
El were longer than those for the channel source used in E2 because the waiting room had 
less surface area to cool the helium/cluster mixture and it emptied slowly through its small 
nozzle. The channel source blew its helium/cluster mixture quickly through its channel 
while providing more rapid cooling. Our final moment measurements were all conducted on 
clusters with residence times sufficient to ensure that they had reached thermal equilibrium 
with their sources and had vibrational temperatures Tyib equal to their source temperatures. 

B. Measured Magnetic Moments 



In both experiments, El and E2, we obtained deflection profiles for each cluster size at 
many temperatures Tvib and many magnetic fields B. We calculated a /iexpt from each profile 
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Number of cobalt atoms 

FIG. 2: Magnetic moments per atom fi for cobalt clusters consisting of 13-200 atoms. Values 
obtained in Experiment El are represented by open circles and those obtained in Experiment E2 
are represented by shaded triangles. The dashed line indicates the single average value measured 
for Co65-.2i5 by Douglass, et al. (Ref. Iioi). 

and then fit all the measurements for a given cluster size to Eq. ([T]). The high quality of 
the fits demonstrates the validity of the Langevin magnetization relationship. The only free 
parameter in each fit is fi, the true magnetic moment per atom for that cluster size. 

Figure [2] shows the values of fi obtained for each cluster size in the two experiments, 
El and E2, and those values are also listed in Table [U The listed values include only the 
statistical uncertainties, which are small and demonstrate the excellent agreement between 
Eq. ([T]) and the experimental observations. The systematic uncertainties are larger and 
reflect limitations in the magnet calibration, position control of the cluster beam in the 
magnet, and the control of vibrational temperature. We estimate this systematic uncertainty 
at ±7% in El and ±5% in E2. 

IV. DISCUSSION 

In addition to the two experiments reported here, there have been two other experimental 
studies of magnetism in free cobalt clusters published recently: a study by Xu et al— 
(hereafter referred to as Xu) covering the entire range of cluster sizes discussed in the present 
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FIG. 3: Magnetic moments per atom fi for cobalt clusters consisting of 13-32 atoms. Points 
marked by hollow circles and shaded triangles correspond to the results of Experiments El and 
E2, respectively. Points marked by shaded diamonds are the results of Xu et al. from Ref. 12 and 
points marked by shaded squares are the results of Knickelbein from Ref. y. 

work and a study by Knickelbein^^ focusing on clusters consisting of between 7 and 32 
atoms. There are many important similarities and differences between the results of the 
four experiments. 

One of their most notable and physically important similarities is that they all find en- 
hanced magnetism, relative to that of bulk cobalt (1.7 yU^/atom), in all of the cobalt clusters 
that deflected toward strong field. This enhancement has been observed before^i^i^i^i^ and 
anticipated theoretically due to the reduced dimensionality of the clusters.""^ Their small size 
leads to an increase in the density of states at the Fermi level and to increased contributions 
from unquenched orbital and spin magnetism.— 

Perhaps the most striking difference between the studies appears in the smallest sizes 
observed. Figure [3] shows values of /i for clusters consisting of between 13 and 32 atoms for 
all four experiments. While our experiments El and E2 found similar /i values, they differ 
significantly from those obtained by Knickelbein and by Xu. Not only do few of the /i values 
agree between these three groups, there is also little correlation between relative maxima 
and minima as a function of cluster size. And while our experiments found /i decreasing 
with increasing size up to C026, Knickelbein found approximately constant up to that size 
and Xu found /i increasing with increasing size over that range. 

While these difference could simply reflect problems with one or more of the experi- 
ments, they could also stem from more exciting thermal and isomeric issues. Both Xu and 



10 



Knickelbein studied clusters smaller and colder than we studied, and thereby observed some 
interesting behaviors. Knickelbein measured a marked change in the deflection profiles of the 
smallest clusters at low temperatures. These clusters apparently shift from a superparam- 
agnetic state at higher temperature to a non-superparamagnetic state at lower temperature. 
This is not surprising because studies of supported clusters have shown that the magnetic 
anisotropy energy per cluster increases dramatically with decreasing cluster size for clusters 
of fewer than 25 atoms.— 

Both Xu and Knickelbein have also reported the presence of two distinct populations 
of magnetic isomers in their beams under certain conditions.—"^ Although geometric iso- 
mers of cobalt clusters have been predicted and studied theoretically,— i^i^i^ the calculated 
differences in their magnetic moments were not as large as those observed in these two 
experiments. 

We found no evidence of a second, less-magnetic population in either of our experiments, 
so the isomer populations in our cluster beams may have been different from those in the 
beams of Knickelbein and Xu. When smaller moment isomers are present in a beam and not 
resolved from the higher moment isomers, they will reduce the average magnetic moments 
and lead to lower reported values of fi. If unresolved low-moment isomers are present in the 
experiments of Knickelbein and Xu but absent in our experiments, that could explain why 
the values for fi that Knickelbein and Xu report are smaller than those we measure. 

For cobalt clusters with up to 26 atoms, theoretical studies^*^!^ have predicted that fi 
should decrease with increasing cluster size, consistent with the trend we observe in El and 
E2. One of those theoretical studies^ also predicts a sharp minimum in fi for C026, a feature 
that appears in both of our experimental results. However, the average /i values predicted 
by these calculations are lower than our average /i value and closer to those of Xu and 
Knickelbein. Nonetheless, the models used in those calculations don't consider unquenched 
orbital magnetism, and recent work has shown that electron orbital motion can contribute 
substantially to the overall magnetic moments of small clusters.— 

For clusters consisting of between 26 and 32 atoms, the /i values reported by Xu are similar 
to our yU values, while those obtained by Knickelbein are still significantly lower. Knickelbein 
didn't study clusters consisting of more than 32 atoms. All three of the experiments that 
continue beyond C032 observe a local maximum in fi at C036-C037, followed by a rapid 
decrease in /x as the next 4 or 5 atoms are added. That sudden dechne in fi suggests that a 
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basic structural change — probably a rearrangement of the underlying geometrical lattice — is 
occurring in the cobalt clusters over this size range. Yang et ai- and Parks et alM also 
found evidence for a structural change in this size range in their ionization potential and 
chemical reactivity studies, respectively. Another study done recently in our laboratory 
suggests that free chromium clusters undergo a striking geometrical change at a slightly 
smaller size (Cr32 to Cr34).— While the theoretical studies of cobalt clusters do not predict 
a decline in /i, they do not consider every cluster size for these larger clusters.— 

Above about C042, all three experiments, and the earlier one of Billas, et al.,— observed 
small oscillations in /i with increasing cluster size (Fig. HI). This behavior is typical of clusters 
of ferromagnetic materials and has been explained in terms of the shell model of cluster 
growth.— In the shell model, the spectrum of cluster sizes is punctuated occasionally 
by unusually stable "closed shells," where what has closed is either a geometrical layer of 
atoms to form complete facets or a set of electronic levels. Long associated with "magic 
numbers" in cluster population distributions, closed-shell clusters are also expected to be 
local minima for fi in ferromagnetic clusters. Beyond each shell closing, the value of fi is 
expected to drift upward initially before decreasing toward a new minimum as the next larger 
shell closes. That behavior was proposed more than a decade ago,— i^i^i^ but extensive 
calculations on iron clusters have recently bolstered the notion that these variations in n 
are due to geometrical factors.^ In their calculations of different lattice structures, Tiago, 
et al. showed that faceted arrangements correspond to local minima in magnetic moment 
per atom fi. 

The differences in the oscillations observed in the three experiments (El, E2, and Xu) 
illustrate how difficult it is to make ideal measurements. Not only are there statistical and 
systematic uncertainties in the measurements, there are also inhomogeneities and irrepro- 
ducibilities in the cluster samples themselves. In each experiment, the cluster samples must 
inevitably include minor isomers: clusters with equal numbers of atoms but slight differences 
in how some of those atoms are arranged on the cluster surfaces. There is now evidence 
that more significant isomers are also present,—*^ isomers with very different magnetic mo- 
ments and properties and with relative populations that may depend on temperature and/or 
carrier gas pressure. 

One noticeable difference between our experiments and that of Xu is that our /i values 
are consistently about 0.2 /i^ per atom larger than Xu's. This difference could be due 
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FIG. 4: A review of the measured magnetic moments per atom /i for cobalt clusters, including both 
the present work and all of the previous measurements of beam-isolated cobalt clusters. Points 
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to systematic errors in the measurements themselves. First, Xu's clusters could be hotter 
than expected, or ours colder. Because each experiment includes data collected at several 
different temperatures, however, one or both experiments would have to be mistaken about 
temperature by the same percentage at each temperature studied. 

Second, a systematic error in either lab in the calibration or alignment of the gradient- 
field magnet, or in the measurement of the beam position in that magnet, could explain 
the discrepancy. To reduce the possibility of such a systemic error in our lab, we carefully 
reexamined the beam position and alignment between experiments El and E2, using a 
completely different approach from the one used before El. It is thus unlikely that errors in 
our magnet position or alignment are responsible for the discrepancy, although inaccuracy 
in our magnet calibration remains a possibility. 

There is also an interesting trend to the small differences between the /i values measured 
in our two experiments. El and E2. For C030 to C065, the measured values of /i were typically 
a little lower in El than in E2, while for Coes up, the measured values were typically a little 
higher in El than in E2. This switch could be the result of a gas expansion effect in El. In 
experiment El's "waiting room" source, gas leaving the source at long residence times may 
have cooled itself by doing work on gas leaving the source at short residence times. This 
same cooling effect occurs when you open a can of compressed air. 

Since larger clusters take longer to form and cool, we usually study them at longer 
residence times than the smaller clusters. In El, those larger clusters may have equilibrated 
with gas that had cooled itself slightly below the actual temperature of the source and those 
larger clusters may be slightly colder than we expect. Overestimating their temperatures 
when fitting data from these clusters to Eq. ([1]) would lead us to overestimate their magnetic 
moments per atom. In E2, the helium's ongoing thermal equilibration with the long narrow 
channel through which it flowed should have prevented this sort of excess cooling. 

We looked for a cooling effect in the El residence time studies themselves, but could 
not prove its presence or absence convincingly. While the residence time study in Fig. [U^a) 
may include a slight peak in /i near 1400 fis, consistent with a gas temperature that has 
dropped slightly below the source temperature, the experimental uncertainty is too large for 
any definitive observations. Still, it would not take much excess cooling in El to explain the 
small differences in fi between El and E2 for Cogs up. It is worth noting that almost all of 
the fi values obtained in El and E2 for Cogs up fall within the 2.24±0.14 /i^/atom range 
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reported a decade ago in Ref. Il9| for Coes through C0215. 

We have interpreted the data from experiments El and E2 using Eq. ([1]), which we 
obtained from the superparamagnetic modeL However, the apphcabihty of that model 
to small, isolated clusters remains a matter of some controversy. In its simplest form, the 
superparamagnetic model predicts that clusters of a given size should deflect homogeneously 
to strong field so that there should be no broadening of the deflected beam profile. After 
observing substantial broadening in the deflected profiles, Xu et al. proposed an alternative 
model for the magnetic response of the clusters. In their adiabatic magnetization model, 
clusters entering a magnetic field shift adiabatically into magnetic states through avoided 
crossings between Zeeman levels.— i^i^ 

According to this adiabatic magnetization model, clusters deflect towards strong field on 
average, but diversity of quantum states in the ensemble of clusters entering the magnetic 
field causes the deflected beam profile to broaden significantly. Despite that broadening, Xu 
et al. report that the average magnetization follows the Langevin function and Eq. ([1]),— so 
that the adiabatic model produces average deflections that are indistinguishable from those 
produced by the superparamagnetic model. 

Xu et al. contend that the presence of broadening in the deflected beam profile is incon- 
sistent with the superparamagnetic model because the only explanation for such broadening 
in superparamagnetic clusters is imperfect time averaging — the clusters must fail to explore 
the entire Boltzmann distribution during the measurement time in the magnet, tmag- When 
the thermal relaxation time r is short compared to tmag, time averaging is complete and 
superparamagnetic theory predicts that equal clusters will exhibit equal magnetizations. 
A beam of such clusters will deflect without broadening. But as r approaches tmag, time 
averaging becomes imperfect and equal clusters will begin to deflect unequally. 

For incomplete time averaging to be responsible for the broadening observed in the de- 
flection profiles, r must be both remarkably long and fortuitously close to tmag- Super- 
paramagnetic theory predicts that the width of the magnetic distribution AM/ /i should be 
approximately equal to (T/tmag)^^^- Xu et al. found /S.M/ jj, to be insensitive to tmag, as well 
as to T and to N, and therefore rejected superparamagnetism as the explanation for the 
Langevin-like magnetization behavior. 

However, imperfect time averaging is not the only explanation for profile broadening in 
superparamagnetic clusters. In fact, imperfect time average is the least likely cause of such 
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broadening, r is almost certainly in the subnanosecond range,-i^i2i so that (T/tmag)^''^ ~ 10~^ 
for most cluster beam deflection experiments. Superparamagnetic clusters thus explore the 
Boltzmann distribution so quickly compared to the measurement time that time averaging 
is essentially perfect and equal clusters exhibit equal magnetization. That equality, however, 
does not rule out broadening caused by the spatial structure of the deflecting magnetic fleld 
or by true inequalities between clusters. 

The spatial structure of the deflecting magnetic field is clearly responsible for some of the 
broadening. Equation ([T]) shows that fies, the time-average projection of the cluster's true 
magnetic moment onto the magnetic fleld, is approximately proportional to the magnitude 
of that magnetic fleld. Since the cluster beam has a flnite width and the magnet through 
which it passes has a gradient magnetic fleld, clusters on the low fleld side of the beam will 
have smaller /igg and deflect less than clusters on the high fleld side. The magnetic fleld 
gradient also varies along the height of the cluster beam, perpendicular to the plane of the 
experiment, further broadening the experimental deflection proflles. 

The second and most important cause of broadening are the inequalities between clusters 
resulting from the thermal statistics of small systems. The fi^s given by Eq. ([1]) depends on 
temperature, but what is the temperature of a particle with only a few dozen atoms? This 
question is especially appropriate in light of the fact that isolated clusters have several dif- 
ferent temperatures (translational, vibrational, and rotational) that are often not in thermal 
equilibrium with one another.— 

Thermodynamics is based on the fact that at a temperature, T, a system's extensive 
properties will fluctuate with a variance on the order of where N is the number of 
particles in the system. It was designed to describe phenomena in macroscopic systems, in 
which is very large and the energies associated with those phenomena are much greater 
than ksT. In nearly inflnite systems, fluctuations produce insigniflcant and essentially 
undetectable fluctuations in the extensive properties. 

But fluctuations are not insigniflcant in flnite systems such as clusters. Recent work has 
shown that if the energies associated with phenomena in flnite systems are of the order 
ksT, then the fluctuations in the value of an extensive property can be on the order of that 
property's average value.— "^"^ 

In our experiment, the clusters' magnetic energies NfiB are rarely larger than fcsT. For 
example, for Coioo at 100 K and B = 0.15 T, values typical of this experiment, NfiB ~ 
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0.2A;bT. Fluctuations in observed magnetism should therefore be substantial and could be 
on the order of those experimental moments themselves. Such statistical distributions in 
experimental magnetic moments, together with the likely presence of minor isomers, should 
cause the width of a beam of superparamagnetic clusters to increase significantly as it 
deflects. 

We studied deflection broadening as a function of all the experimental parameters and 
found empirically that Arf, the broadening of a deflection profile's full width at half maxi- 
mum, is proportion to d, that cluster's average deflection, divided by the number of atoms 
in the cluster (Fig. [5]). That relationship can be written as 

- - - (2) 

d ~ AT' 

where C is a constant for cobalt clusters. The slope of the line in Fig. [5] yields C ~ 54. 

In Stern-Gerlach cluster deflection experiments, the experimental moment per atom /igxp 
is related to the deflection d according to 

Tndv'^ 

where m is the mass per atom, v is the cluster beam velocity, V-B is the field gradient, L is 
the magnet length, and D is the drift region length. Substituting Eq. into Eq. ([2]) gives 

A/^exp _ C 



/^exp ^ 

Far from saturation, Eq. ([T]) can be approximated as 



(4) 



fieS ^ ■ (5) 



Assuming that fies in Eq. ([5]) is /icxp, Eq- © becomes 

AT^ib C 



(6) 



Tvib iV' 

where ATvib is defined as a positive quantity in order to eliminate a meaningless negative 
sign. 

This analysis suggests that the clusters have a fractional spread in vibrational temperature 
ATvib/Tvib that is proportional to C/N. That this fractional spread increases with decreasing 
cluster size as is consistent with the statistics of finite systems and the inevitable 
deviations from the equipartition theorem. 
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FIG. 5: As cobalt clusters consisting of N atoms deflect toward strong field in our apparatus, 
their deflection profile broadens and its full width at half maximum increases. Plotted here are 
experimentally measured broadenings in profile width (i.e., the difference between field-on width 
and field-off width) as a function of average field-on deflection divided by N. The line is a least- 
squares fit to the experimental points and has an approximate slope of 54. 

One hallmark of superparamagnetic behavior is that the clusters should never deflect 
toward weak field, even when broadening effects are taken into account. The adiabatic 
magnetization model, on the other hand, does not exclude weak-field deflection. We looked 
carefully for deflection to weak field but could not find it in any of our measurements. 

In fact, when we studied the deflection of the cluster beam's low-field edge by strong 
magnetic fields and field gradients, we found that we could deflect that edge of the beam 
toward strong-field by as much as 1 mm. This deflection was complete in that there were 
no clusters left undeflected at the beam's low-field edge. If any of the clusters had deflected 
toward weak-field or avoided deflection altogether, there would not have been an empty 
region along the low-field side of the beam. We found no clusters there, indicating that all 
of the clusters deflected toward strong-field, consistent with superparamagnetism. 

Although both cluster magnetism models, superparamagnetism and adiabatic magneti- 
zation, yield Langevin-like relationships between the effective moment per atom /igg and 
the true moment per atom /i, the temperature that is relevant to the Langevin function is 
different for the two models. Clusters in a seeded supersonic beam have at least three dif- 
ferent temperatures, a translational temperature Ttrans, a rotational temperature Trot, and a 
vibrational temperature T^ib, and those three temperatures are unlikely to be in equilibrium 
with one another. 
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It is well established that when molecules or clusters are seeded in a noble gas and undergo 
supersonic or free jet expansion, their translational and rotational temperatures cool far 
more easily than their vibrational temperatures.— i^i^i^ In most cases, Tt^ans < ^rot < ^vib- 
For example, Wallraff et al. measured Tj-ot = 5.7 K and T^ib = 213 K in a beam of the 
triatomic molecule OCS from a free jet expansion in room-temperature argon. '^^ While it is 
relatively straightforward to produce a supersonic beam of clusters with Tvib — Tsource,^° it is 
much harder to achieve Trot — ^source- Rotational cooling is so efficient in seeded supersonic 
expansions that clusters in the resulting beam will almost certainly have Trot << ^source- 

A second important thermal issue is the heat capacity associated with each cluster tem- 
perature. Xu et al. assert that a cluster's vibrational degrees of freedom should freeze out 
when Tvib < ^Dcbye/A^^^^ and therefore that most of the smaller clusters they study are in 
their ground vibrational states. This suggestion that the clusters effectively have zero vi- 
brational heat capacity at the experimental temperatures neglects the softening that occurs 
when a solid is reduced to cluster dimensions. With most of its atoms on its surface, a 
cluster is no longer bulklike with respect to phonons and vibrations,—*^ as evidenced by the 
reduced melting temperatures-^^ and premelting effects^ seen in nanoparticles and clusters. 
Instead, it retains considerable vibrational heat capacity even at relatively low temperatures. 
For example, carbon has a Debye temperature of 2230 K, yet Ce has vibrational energies as 
low as ~ 100 cm-i (- 140 K).^^ 

In superparamagnetism, each cluster acts as its own heat bath and that heat bath is 
dominated by the cluster's vast number of vibrational degrees of freedom. Moreover, the 
cluster's three rotational degrees of freedom are highly constrained by angular momentum 
conservation. As a result, the effective temperature of the cluster's internal heat bath and 
the temperature that must appear in the Langevin function is Tvib. Fortunately, Tvib is 
relatively easy to control in a seeded supersonic cluster beam. 

In the adiabatic magnetization model, however, the assumption that spin-rotation cou- 
pling is responsible for the Langevin-like behavior implies that the rotational temperature 
Trot is dominant. Moreover, if the vibrations are frozen out, Tvib has no role at all in this 
model. The temperature that should appear in the Langevin function is Trot- Since Trot is 
difficult to control, and unlikely to be even close to the Tgource, measurements of magnetic 
properties that depend on Trot will be difficult to interpret. 

Each cluster has a geometrical lattice and the cluster's magnetic moment tends to align 
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along an easy axis of that lattice. The energy required to reorient the moment away from an 
easy lattice axis is known as the magnetocrystalline anisotropy energy (MAE). A familiar 
concept for bulk materials, MAEs have also been calculated for cobalt clusters and are found 
to depend strongly on the size and detailed structure of a cluster.— For cobalt clusters 
of the sizes we studied, the MAEs are predicted to range from approximately 30 /ieV (0.4 
K) for highly symmetric and/or small clusters to 400 /ieV (5 K) for less symmetric and/or 
large clusters. Since Tvib > 60 K for all of our measurements, ksT is large compared to the 
MAE and our clusters should be well into the superparamagnetic regime. 

For cobalt clusters to avoid superparamagnetic fluctuations, they would have to be cooled 
to extremely low temperatures or have anomalously large MAEs. Knickelbein reports ob- 
serving non-superparamagnetic or "locked-moment behavior"^- in Co7_ii at temperatures 
below 100 K and in C015 at temperatures near 50 K.— 
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s also have a role in the adiabatic magnetization model proposed by Xu, et al. in 
In that model, the Zeeman levels arising from a cluster's magnetic spin quantum 
number 5* fan out in energy with increasing magnetic field. Since each rotational level 
associated with the cluster's rotational quantum number J has its own manifold of Zeeman 
sublevels, the Zeeman fans of different rotational levels encounter one another as the field 
increases. Couplings between spin and rotation cause these Zeeman levels to experience 
avoided crossing. The adiabatic magnetization model proposes that as each cluster enters 
a magnetic field, its state evolves adiabatically through the avoided crossings in such a way 
that an ensemble of clusters exhibits an average magnetization that depends on the Langevin 
function Eq ([T]). 

The spin-rotation coupling responsible for these avoided crossings is the MAE and without 
that coupling, the Zeeman levels would simply cross one another. But this Zeeman-level 
picture assumes that spin-rotation coupling is relatively weak compared to the level spacings 
so that a cluster's spin angular momentum 5* and rotational angular momentum J are both 
good quantum numbers. 

In most of the experimental conditions, however, the rotational energy spacings are less 
than the MAE. For example, at 60 K, C055 has J ~ 270 and the energy spacing between J 
and J + 1 is approximately 30 /ieV. Since the predicted MAE for C055 is also approximately 
30 /ieV, the coupling energy is comparable to the spacings between rotational levels and the 
Zeeman picture blurs into a sea of mixed states. For larger or less-symmetric cobalt clusters. 
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the spin-rotation couplings are even larger. A detailed theory of these spin-rotation effects 
should therefore be rather complicated. 

V. CONCLUSION 

We have measured the magnetic moments per atom for small cobalt clusters (13-200 
atoms) and find them size dependent. All of these clusters have magnetic moments per 
atom that are larger than the value for bulk cobalt (1.7 /i^/atom), evidence that the re- 
duced dimensionality and increased surface to volume ratio of clusters leads to enhanced 
magnetism. We see gradual magnetic oscillations with increasing size, consistent with a 
shell model of cluster growth. Changes made to our experimental apparatus and data 
acquisition procedure were found to have had little impact on the results. We find that su- 
perparamagnetism and the statistics of finite systems can explain both our experimentally 
observed magnetic moments per atom and the spreads in the deflection profiles of all of the 
clusters. We have compared the superparamagnetic and adiabatic magnetization models 
with the experimental and theoretical evidence and find that the superparamagnetic model 
is the most likely explanation for the Langevin-like magnetizations of the isolated cobalt 
clusters in our experiments. 
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TABLE I: Magnetic moments per atom, /x, obtained by ap- 
plying Eq. (pP) to the experimentally measured moments 
per atom, /Xeff. Only statistical uncertainties are indicated 
(1 SD). 
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